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Tensor networks representations of many-body quantum systems can be described in terms of 
quantum channels. We focus on channels associated with the Multi-scale Entanglement Renormal- 
ization Ansatz (MERA) tensor network that has been recently introduced to efficiently describe 
critical systems. Our approach allows us to compute the MERA correspondent to the thermody- 
namic limit of a critical system introducing a transfer matrix formalism, and to relate the system 
critical exponents to the convergence rates of the associated channels. 
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Understanding the properties of strongly interacting 
many-body quantum system is central in many areas of 
physics. Whenever it is hard to device rehable analyt- 
ical approaches, as in many situations of experimental 
relevance, ingenious numerical methods are necessary to 
grasp the essential properties of these systems. Our abil- 
ity of simulating them is based on the possibility to find 
an efficient description of their ground state. This is the 
case, for example, of White's Density Matrix Renormal- 
ization Group [l| which can be recasted in terms of Ma- 
trix Product States (MPS) 0, i, i, i, i]. Such repre- 
sentations are characterized by a simple tensor decom- 
position of the many-body wave-function which allows 
one i) to efficiently compute all the relevant observables 
of the system (e.g. energy, local observables, and cor- 
relation functions), and ii) to reduce the effective num- 
ber of parameters over which the numerical optimization 
needs to be performed. MPS fulfill these requirements 
and can be used to describe faithfully the ground states 
of not critical, short range one-dimensional many-body 
Hamiltonians at zero-temperature. However MPS typ- 
ically fail to provide an accurate description in other 
relevant situations, i.e. when the system is critical, in 
higher physical dimensions or if the model possesses long- 
range couplings. Several proposals have been put for- 
ward to overcome this problem. Projected Entangled 
Pair States (PEPS) 01 generalize MPS in dimensions 
higher than one. Weighted graph states [8| can deal with 
long-range correlations. In this Letter we focus on a so- 
lution recently proposed by Vidal [9[ who introduced a 
tensor structure based on the so called Multiscalc Entan- 
glement Renormalization Ansatz (MERA). The MERA 
tensor network satisfies both the constraints i) and ii) 
and accommodates the scale invariance typical of crit- 
ical systems 
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The relevance of this approach 
might represent a major breakthrough in our simulation 
capabilities Il2l and motivates an intensive study of the 
MERA [11, [3. 

Here we point out a previously unnoticed connection 
between the MERA and the theory of completely posi- 
tive quantum maps [l5[ establishing a link between two 



important areas of quantum information science. This 
allows us to introduce a transfer matrix formalism in the 
same spirit as it has been done for MPS [3|, [Sj, providing 
new tools to compute physical observables using MERA. 
The main outcomes of our work are i) a method for de- 
termining the properties of critical many-body systems 
in the thermodynamic limit and ii) a connection between 
the critical exponents governing the decay of correlation 
functions and the eigenvalues of the MERA transfer ma- 
trix. As a consequence this yields a full characterization 
of the asymptotic properties of one-dimensional critical 
systems. The paper is organized as follows: after a brief 
review of the MERA, we show how the local expectation 
values and correlations functions can be casted in terms 
of concatenated quantum channels. Usinggeneral prop- 
erties of mixing quantum channels [16l.ll7l.ll8l.ll9l| we then 
provide a way for expressing the thermodynamic limit of 
these quantities. 

The MERA tensor network:- Consider a many-body 
quantum system composed by iV = 2" sites of dimen- 
sion d (qudits). Any pure state can be expressed as 

\^) = E ^i/2,---,«ivlCfi,Cf2r-- :C«iv): where for j e 
{I,--- ,7V} and i e {!,•■• ,d} the vectors \^i)j £ Hd 
form the basis of the j-th d-dimensional system compo- 
nent. The MERA representation [91] of \ijj) assumes a 
specific tensor decomposition of T described in Fig. [TJ 
Here the links emerging from the lowest part of the 
graph represent the N physical indices of T. The nodes 
of the graph instead represent tensors. They are di- 
vided in three groups: the type-T ^ ) disentangler ten- 
sors X of elements Xe^f^ represented by the red Xs; 
the type-(' I ) tensors A of elements A"^^^^ represented by 
the blue inverted Ys; and the type-T ° ) tensor C of el- 
ements Cg^^i^^i^^e^, represented by the green semi-circle. 
As shown in Fig.[T]the x's, the A's are coupled together to 
form a triangular structure with C as the closing element 
of the top: any two joined legs from any two distinct 
nodes indicate saturation of the associated indices. Con- 
sequently, the tensor T associated with the A^ qudit state 
\i/j) is written as a network of 0{N) tensors organized in 
0(log2 N) different levels composed by one layer of x ten- 




FIG. 1: Representation of a MERA decomposition of the ten- 
sor T associated with a many-body state \xp) for A*' = 16 qu- 
dits (a dotted link emerging from the left side of the graph 
re-enter the figure as the corresponding dotted hnk on the 
right). The light blue region represents the causal cone Q] 
associated with the local operator Oj (black rectangle). It is 
composed by contraction of tensors AI5 (shown in the inset). 



sors connected with one layer of A tensors. In a generic 
MERA the x's and the A's differ from node to node but 
the dimensions of their indices are upper-bounded by a 
fixed constant (for easy of notation we omit the labels 
expressing the position of the tensors in the network). 

What makes the MERA decomposition a convenient 
one is the assumption that the x's and A's satisfy spe- 
cial contraction rules. Specifically, for each x f^nd A let 
us define its adjoint x ^^^ ^ ^^ the tensors of elements 
{X":;r:r and A;^; '"^ = (Af/^ ,„ J* . with this def- 

i,-ui^^2,-u2 and A*;* X^^,, = 81^,^, with 
5 being the Kronecker delta and where the typographic 
symbols • and <> indicate summation over the correspond- 
ing index. Under these conditions the expectation val- 
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ues of local observables on |-0) requires only to evaluate 
0(log2 N) non trivial tensor contractions [91]. 

Local observables and quantum channels:- We first 
show how the average of local observables can be related 
to the study of concatenated quantum channels. Given 
an operator 0^ which acts not trivially on no more than 
three consecutive qudits (say the {j — l)th, jth and 
(j -I- l)th), the quantity (Qj) = {ip\Qj\ip) requires to 
perform contractions only over the x's and A's belonging 
to the causal cone [9| of the triple i — 1, j and j -I- 1. 
A compact expression is obtained by grouping these 
tensors in compounds composed by 2 x's and by 3 A's 
(see inset of Fig. [1]). This forms m = log2(A^/4) non nec- 
essarily identical type-T ^ ) tensors M^ = Ax Ax A 
where the products Ax and X'^ are defined by 
[^X£:l\ - A-. x;a, and [xA]--,3 ^ xTJ. K^,- 
For each one of the m tensors A^s we can then in- 
troduce two families of operators {£r}r and {Rr}r 
acting on the Hilbert space Ti.®^ and labeled through 



the composed index r = {ri,r2,r^) with ri2,3 being 
d-dimensional [20|. In the computational basis they are 
defined by the matrices (C«i,?«2i'C«3l-^r-|Cfi,C£2'63) 
and {£.m,£.u2,iu3\Rr\S.f.i,£.i2,£,i3) of elements 

[■^5];\t^:%..2..3 and [M,Trl^;;:;-X,i.,r, respectively. 
They are related through a reshuffling of the input and 
output qudits, i.e. Rr = n(L,.) = PLrP\ where P = P^ 
is the unitary transformation which exchanges the first 
and the third qudit. Most importantly, according to 
the contraction rules defined previously, they satisfy the 
normalization conditions J2r ^^LI = /®^ = ^^RrRl, 
with / being the identity operator of Hd- This implies 
that {Lr}r can be used to define a completely positive, 
unital, not necessarily trace preserving super-operators 
^H uM^ which transforms the linear operator O of 
Hf^ into $^^0) = Er^rQLl- Analogously {Rr}r 
defines the map $)j which is related with ^\^' through 
the identity $^ -* = 11 o <1>^^ o H, where "o" indicates 
the composition of super-operators. We also introduce 
the vector of Hf^, \hat) = EC^i.^jAA \^ii,£.e2,^i3,£.ei), 
which without loss of generality is assumed to be nor- 
malized, and define pc the three sites reduced density 
matrix obtained by tracing \hat){hat\ over one of the 4 
qudits. With these definitions one can finally write the 
expectation value of Qj as (8j) = Tr[/5c -B, ], with 
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<f>)^ (Oj), and where (enumerating 
from the lower MERA level of Fig. [J) <^P is either 

the map <I>)^ or $}j associated with the fc-th tensor 
A^5 of the causal cone (which one depends upon N and 
j). The operator B^^ is thus obtained by applying 

to the observable Qj a sequence of m super-operators 
associated to the MERA causal cone. We can then write 



(e,) = Tr[<i> 
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where $(1*-™) = $(1) o • • • o $(™), with ^C^) being the 

(k) 

super-operator <I>)^ in Schrodinger picture. By con- 
struction the $^*''' (and hence $(1^™)) are Completely 
Positive, Trace Preserving (CPT) maps, i.e. quantum 
channels with Kraus operators [15[ defined by either the 
set {Ll}r or {Rl}r- Equation ([T]) establishes a formal 
equivalence between the MERA tensor network and the 
successive application of a family of CPT maps (the 
QuMERA family). Since it holds for all the local ob- 
servable Qj this implies that ^'^^^'"^(pc) coincides with 
the reduced density matrix pj of the input state \ip) as- 
sociated with the qudits j — 1, j and j + 1, i.e. 



$('^'"^(/5c) = 



Pi 



(2) 



Critical systems in the thermodynamic limit: - Con- 
sider now a family of MERA states ji/j) which describes 
the ground state of a many-body Hamiltonian H at 
criticality. We are specifically interested in the ther- 
modynamic limit of infinitely many sites (i.e. N — > 



cxd). According to the above derivation, in the hmit 
of large m Eq. ^ converges toward the reduced den- 
sity matrix px of three consecutive qubits of the systems 
hmm^oo *&^^^'"'(Pc) = Pt- Of course the above Hmit 
should not depend upon the particular causal cone "tra- 
jectory" one chooses to follow (the system is translational 
invariant). Without loss of generality we can thus pick 
the one associated with the central sites of the MERA, 
i.e. the one associated with the causal cone of 7V/2-th 
qudit. This allows us to identify aU the ^^''^ of $(1^™) 
with maps of the form $(^). A further simplification 
arises by enforcing the scale invariance property of the 
system. This can be done for instance by assuming that 
all the tensors x's and A's of the MERA to be identi- 
cal III and by requiring $(^) = $(^) = $ [2l[. With this 
assumption all the sequence $(i^™) can now be written 
as a composition of m identical quantum channels, i.e. 
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FIG. 2: Causal cone associated with the two point correlation 
function (Oi©-,). The light blue regions generate the tensor 
product channel ^t^^™-' (g $(1^™). The magenta region in- 
stead corresponds to the CPT map \]/('"*~'") which, acting on 
\hat) produces the 6 qudits state p\"^ ■ 



By general results on quantum channels the vast majority 
of CPT maps are known to be mixing (or relaxing) |16l.ll7l. 
|l8|, [l9|. This means that for a generic choice of $, in the 
limit m —>■ 00 the transformation ([3]) will send all input 
states into a unique fix point identified as the unique 
eigen-operator of $ associated with its largest eigenvalue. 
This property allows us to identify the thermodynamic 
limit Pt of the reduced density matrix 5_; with such an 
eigenstate. As in the case of MPS [3, [23, we can now 
provide a simplified expression for the thermodynamic 
limit of any local observable & of scale invariant MERAs. 
A convenient way to express this is obtained by moving 
in Liouvillc space |15l . |23[ . By doing so we can write 



{&)t = lim {{Q\{E^r\pc)) = mpr)) 



(4) 



where the vectors \A)) are the Liouvillc representations 
of the operators A, while E^ is the Liouvillc operator 
associated to the map $, i.e. 



i?<J> ^ y Rr ' 



'R: 



(5) 



(here R* is the complex conjugate of Rr evaluated with 
respect to the canonical basis). Within this formalism the 
vector \pt)) which describes the fix point pT of the mix- 
ing map $ is also found as the the (unique) eigenvector 
correspondent to the unitary eigenvalue of E^ . Because 
of the close similarities between E^ and the MPS trans- 
fer matrix [2, Isj, la, 122] we dubbed the former the transfer 
operator of the MERA. The existence and uniqueness of 
the fix point is given by the physical assumption that the 
thermodynamical limit of a physical system exists and it 
is unique. 

Correlation functions:- The computation of the long 
range correlation functions of {ip) has also a clear inter- 
pretation in terms of CPT maps and the thermodynamic 



limit can be computed along the same lines presented for 
local observables. Here we specialize in the two point 
correlation functions as the generalization is straightfor- 
ward. Consider then the expectation value (QiQj) with 
Qi and Qj being two local operators acting on (say) the 
i-th and j-th qudit respectively. In this case the causal 
cone is formed by two single sites causal cones which in- 
tercept at the MERA level fh = int[log2(J — j)] — 1, (see 
Fig. 12]). The resulting expectation values can then be 
written as 



(9,9,) = Tr[($ 
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(1- 



')(p!f)(e,0 9,)], (6) 



where $,- *~"^' are the two CPT maps of the two single- 
site causal cones associated with the sites i and j re- 
spectively (light blue regions of Fig. l2]). The operator 
Pj-™ instead is a 6 qudits state associated with the last 
m — fh levels of the MERA. It is obtained from the 4 
qubit state \hat) through the application of a quantum 
channel \I'(""~™) which, similarly to ^l I~"^\ originates 
from a proper concatenation of CPT maps associated 
with A^5 or with the type-( ^^^ ] and type-( ^ ^ tensors 
A4g = AxAxAxAxA and Air = AxAxAxA. Since this 
applies to all the two sites observable, we can then con- 
clude that {^f™' (g) ^y'^"^'){pl™') must coincide with 
the reduced density matrix pij of \hat) associated with 
the sites i and j. 

Let us focus then on the thermodynamic limit of the 
correlation function Ay = (9^9,) — (Qi){Qj) which for 
Hamiltonian systems at criticality decays as \i — j\~'^ ■ 
Under the same assumptions used to derive Eqs. ([3]) and 
(|4|) we can assume j to be the central site of the MERA 
(i.e. j = N/2). Suppose then that the associated map <& 
is mixing with fix point pT- For any input state pij of the 



sites i and j we then have Hnim^oo (2'i ® $^ ^"^ ) (/5y ) = 
Pi (S) PT, with Pi = Tij[pij] and 2^ being the identity 
super-operator of the site i. The speed of convergence, 
evaluated through the trace distance, is exponentially 
fast [l8|, ll9[ in fh and, a part from some constant pre- 
factor, can be upper-bounded by the quantity fh''' k™ 
with K < 1 being the modulus of the largest eigen- 
value of $ whose associated eigenvector contribute in 
the expansion of Eq. ^. This is sufficient for claiming 



that the distance between ($, 
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{ ^rnj^^^^ j^ bounded by oc k 



2" Thus 



we can write log2 (|Ay |) ^ 2mlog2 k + C'(log2 m), which 
through the definition of fh provides a bound for the crit- 
ical exponent v associated to the observable in terms 
of the properties of the map $, i.e. v ^ — 21og2 k. In 
effect one can show that such bound is tight, i.e. 



~21og2 K 



(7) 



This can be seen for instance by expressing Eq. ^ in the 
Liouville space formalism as in Eq. (|3|), and expanding 
the transfer matrices -E^. (8) E<^. in Jordan blocks (the 
calculation is similar to the MPS analysis of Ref. [22j). 
A numerical test of Eq. (0 [li] on a MERA state ap- 
proximating the ground energy of an Ising chain up to 
a 10~^ accuracy yielded k ~ 0.915, 0.49, 0.52 to be com- 
pared with the exact values Kth — 0.917, 0.46, 0.50 associ- 
ated with the X, y and z two-point correlation functions. 
Similar results have been obtained for the XXZ model. 

Concluding remarks:- Equations. (O and ([7]) consti- 
tute the main results of our analysis. As already dis- 
cussed by Vidal [9| , MERA networks are able to describe 
algebraic decaying correlations. In this work we put on 
firm grounds this observation giving an explicit expres- 
sion of the critical exponents in terms of properties of 
the associated QuMERA channels. The combination of 
this approach with conformal field theory methods may 
provide a powerful tool to achieve a complete descrip- 
tion of one-dimensional critical quantum systems. Simi- 
larly our findings yield a natural connection between the 
tensor network description of the thermodynamic limit 
of critical systems and the master equation formalism. 
Combining these results with the algorithms presented 
in 10|, Il3| one can exploit the introduction of the trans- 



fer operator ^ studying directly the infinite size system 
improving simulation efficiency. 

The results presented here can be easily extended in 
several ways. For instance since a binary tree can be 
seen as a MERA with the disentanglers x set to the iden- 
tity, all the arguments presented previously can be easily 
adapted to this case. Similarly also the thermodynam- 
ical limit for MPS [5|] can be described in terms of re- 
peated application of CPT maps (here the operator ([5|) 
reduces to the MPS transfer matrix). More generally our 
approach can be adapted to any tensor network by asso- 
ciating it with a family of CPT transformations which. 



properly concatenated, allows one to compute the local 
observables of the system. In this perspective the quan- 
tum circuit [sl, llJ] associated with the tensor network 
can be seen as a unitary dilation or, Stinespring rep- 
resentation [15[, of the corresponding CPT family. Fi- 
nally a generalization to higher spatial dimensions seems 
straightforward. 

We thank M. Rizzi for useful discussions. This work 
was in part founded by the Quantum Information re- 
search program of Centro "Ennio De Giorgi" of SNS. 
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